On the role of transition region on the Alfven wave phase 
mixing in solar spicules 
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Abstract Alfvcnic waves arc thought to play an im- 
portant role in coronal heating and solar wind acceler- 
ation. Here wc investigate the dissipation of standing 
Alfven waves due to phase mixing at the presence of 
steady flow and sheared magnetic field in the strati- 
fied atmosphere of solar spicules. The transition re- 
gion between chromosphere and corona has also been 
considered. The initial flow is assumed to be directed 
along spicule axis, and the equilibrium magnetic field 
is taken 2-diniensional and divergence-free. It is deter- 
mined that in contrast to propagating Alfven waves, 
standing Alfven waves dissipate in time rather than in 
space. Density gradients and sheared magnetic fields 
can enhance damping due to phase mixing. Damping 
times deduced from our numerical calculations are in 
good agreement with spicule lifetimes. Since spicules 
are short living and transient structures, such a fast 
dissipation mechanism is needed to transport their en- 
ergy to the corona. 

Keywords Sun: spicules ■ Alfven waves: phase mixing 
• Transition region 



1 Introduction 

The heating mechanism of the solar cor ona is one of 
the m ost mysterious issue in astrophysics (jAschwanden 
20041 ). In the corona, the temperature rises to a few mil- 
lion degrees. To maintain such a high temperature in 
the corona in spite of cooling by heat conduction and ra- 
diative losses, a continuous supply of thermal energy is 
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necessary. X-ray observations from space experiments 
(e.g. skylab) have shown that the corona is not uni- 
form and consists of many brig ht loops. Two promis - 
ing models for coronal heating (JErdelvi fc Ballaill2007l ) 
are: heating by small scale flares triggered by magnetic 
reconnection, and heating by the dissipation of Alfven 
waves that propagate in the magnetic flux tubes (jAlfven 
1947tlHollwe3ll973i[l986ilMcKenzie et al.lll995l) . How- 



ever, there are other types of waves such as acous- 
tic, slow-mode and fast-mode waves which are strongly 
damped or reflected at the steep density and tem- 
perature gradients of chromosphere and transition re- 
gion. As an origin of Alfven waves, iKudoh fc Shibata 



( 19991) considered a photospheric random motion prop 
agating along an open magnetic flux tube in the solar 
atmosphere, and performed MHD simulations for so- 
lar spicule formation and the coronal heating. They 
have shown that Alfven w aves transport su f ficient en- 
ergy flux into the corona. iDe Pontieu et al.l (J2007l ) es- 
timated the energy flux carried by transversal oscilla- 
tions g enerated by spicu l es and compared by the re- 
sult of iKudoh fc Shibatal (|l999l ). They indicated that 
the calculated energy flux is enough to heat the quiet 
corona and to accelerate the high-speed solar wind. 
There are many proposed mechanisms for the dissi- 
pation of Alfven waves. Interaction of Magnetohy- 
drodynamic (MHD) waves with inhomogeneous plas- 
mas in which has density gradients in both x- and 
z-directions, can leads to a number of very inter- 
esting physical pheno mena, such as mode c oupling 



( Morivasu et al.l 120041) . resonant absorption (jlonson 



Il978t iRuderman et al.l Il997f). magnetohy drodynamic 
(MHD) turbulence (iMatthaeus et al.lll999l ). and phase 
mixing ( Hevvaerts fc Priestlll983l ). These phenomena 
can affect the MHD wave dynamics and cause enhanced 
dissipation and heating of astrophysical plasmas. The 
idea is that when a density gradient of a medium is 
perpendicular to the magnetic field, the Alfven speed 



is a function of the transverse coordinate. On each 
magnetic field hne, Alfven waves propagate with their 
local Alfven speed. After a certain time and distance, 
as height increases, perturbations by Alfven waves of 
neighboring magnetic field lines become out of phase 
and have different wavelengths, causing large gradients 
in Alfven wavcfront in the direction of the inhomogenc- 
ity, that is the phase mixing. Dissipation then appears 
allowing the energy in the wave to heat the plasma. 

Phase mixing of propagating Alfven waves in a strat- 
ified atmos phere, partic u larly, in solar spicules has been 
studied bv Ebadi et all |2012bl ). They concluded that 
dissipation of propagating Alfven waves can occur in 
space rather than in time. Moreover, they found that 
the calculated da mping times are much longer than 



spicule life times. lEbadi fc Hosseinpouij (|2013l ) have 



studied the phase mixing of Alfven waves in the case of 
shear flows and fields. It is determined that the shear 
flow and field can enhance the da mping of standing 
Alfven waves. ISmith et al.l (J2007I ) have showed that 
the enhanced phase mixing can occur both with den- 
sity gradients and sheared magnetic fields. Spicules 
are one of the most fundamental components of the 
solar chromosphere. They are seen in spectral lines 
at the solar limb at speeds of about 20 — 25 km s~^ 
propagating from photos phere into the ma gnetize d low 
atmosphere of the sun ( Zaaarashvili fc E rdclvi 2 0091 ). 
Their diameter and length varies from spicule to spicule 
having the values from 400 km to 1500 km and from 
5000 km to 9000 km, respectively. The typical life- 
time of them is 5 — 15 min. The typical electron 
density at heights where the spicules arc observed is 
approximately 3.5 x 10^^ — 2 x 10^^ m~^, and their 



temp e ratures are est im ated as 5000 — 8000 K ( Beckerg 
1968t ISterli^ I2OOOI). Ikukhianidze et al.l |2006l) and 



Zaaarashvili et al.l (|2007l ) observed their transverse os- 



cillations with the estimated period of 20 — 55 and 
75 — 110 s by analyzing the height series oi Ha spectra 



in solar limb spicules observed. Recently, lEbadi et al 



( 2012al ) based on Hinode /SOT observations estimated 
the oscillation period of spicule axis around 180 s. They 
concluded that the energy flux stored in spicule axis 
oscillations is of order of coronal energy loss in quiet 
Sun. Despite the large body of theoretical and obser- 
vational works devoted to the spicules, their ejection 
mechanism is not clear yet, and observations with high 
spati al resolutions are needed t o dist inguish their ori- 



gin. lOkamoto fc De Pontieaul (J201l[ ) investigated the 



statistical properties of Alfvenic waves in solar spicules 
based on SOT/Hinode observations. They observed 
a mix of upward and downward propagating as well 
as standing waves. They found the occurrence rate 
of 20% for standing waves. They concluded that in 



some spicules there are waves propagating upward and 
downward t o form a standing wave in the middle of 
the spicule. ICargill et al.l (119971 ) performed the numer- 
ical simulations of the propagation of Alfvenic pulses in 
two dimensional magnetic field geometries. They con- 
cluded that for an Alfvenic pulse the time at which dif- 
ferent parts of the pulse emerge into the corona depends 
on the plasma density and magnetic field properties. 
Moreover, they discussed that this mechanism can in- 
terpr et spicule ejection forced th r ough the transition re- 



gion. iMurawski fc Zaqarashvilil (|2010l ) studied the up 



ward propagation of a velocity pulse launched initially 
below the transition region. The pulse quickly steep- 
ens to a shock. They concluded that such a model may 
explain speed, width, and heights of classical spicules. 
This can be a motivation to study the phase mixing of 
upward propagating Alfven waves. 
In this paper we arc interested to study the effect of 
transition region between chromosphere and corona on 
phase mixing. Section 2 gives the basic equations and 
theoretical model. In section 3 numerical results are 
presented and discussed, and a brief summary is fol- 
lowed in section 4. 



2 Theoretical modeling 

Wc consider effects of the stratification due to gravity 
in 2D x-z plane in the presence of steady flow and shear 
field. The phase mixing and the dissipation of standing 
Alfven waves are studied in a region with nonuniform 
Alfven velocity both along and across the spicule axis. 
Non-ideal MHD equations in the plasma dynamics are 
as follows: 



|.V-„v, 



0, 



(1) 



/0^ + /o(v-V)v = -Vp + pgH (VxB)xB + pi^V^v, 

dt fio 

(2) 



dB 

'dt 



V X (v X B) -f ?/V^B, 



VB = 0, 



P 



pRT 



(3) 
(4) 

(5) 



where v and 77 are constant viscosity and resistivity 
coefficients, /ig is the vacuum permeability, fi is the 



mean molecular weight. The typical values for 77 in the 
solar chromosphere and corona are 8 xlO^T~^'^ and 
]^g9j'-3/2 ^ g-i^ respectively. The value of pi/ for a 
fully ionized H plasma is 2.2 x 10~^^T^/^ kg m~^ s~^ 
(jPriestl 119821 ). We assume that spicules are highly dy- 
namic with speeds that arc significant fractions of the 
Alfvcn speed. Perturbations are assumed to be inde- 
pendent of y, i.e.: 

V = vok + Vy{x,z,t)j, 

B = Bo^{x,z)i + Boz{x,z)k + by{x,z,t)j (6) 

and the equilibrium sheared magne tic field is two 



dimensional and divergence-free as ([Del Zanna et al. 
2005HMurawski fc Zaaarashvilill201ol) : 



Box{x,z) 
Bo2{x,z) 



Bqc '"'''cos[A:;,(a; - a)], 
-Soe"'"''^ sin[fc;,(x - a)] 



(7) 



where a is the spicule radius. Since the equilibrium 
magnetic field is force- free, the pressure gradient is bal- 
anced by the gravity force, which is assumed to be 
g=—gk via this equation: 



-\/po{x, z) + po{x, z)g = 0, 

and the pressure in an equilibrium state is: 

dz' 
Po{x,z) ^poix) exp ' 



(8) 



.. A(z') 

The density profile is written in the following form: 

dz' 



Po[x,z) = exp 

To 2; 



. A(^') 



(10) 



where po{x) is obtained from the Alfven velocity for 
a phase mixed and stratified atmosphere due to grav- 
ity which is assumed to be (jDe Moortel et al.l Il999 : 
Karami fc EbrahimillJOOol) : 



po(a;) = /5o[2 + tanh(a(a; - a))] ^, 

and 

RToiz) 



A(z) 



fj-g 



(11) 



(12) 



where po is the plasma density at z = 6000 km, a con- 
trols the size of inhomogeneity across the magnetic field. 
The temperature profile is taken here as a smoothed 
step function, i.e.: 



To{z) - It, 



l + dt + {l- dt) tanh(^ — -) 



(13) 



here, dt = Tch/Tc which Teh and Tc denote the chromo- 
spheric temperature at its lower part and the coronal 
temperature that is separated from the chromosphere 
by the transition region. Zw = 200 km is the width of 
transition region which is located at the zt = 2000 km 
above the solar surface. We put Teh = 15 x 10"^ K and 
Tc = 3 X 10^ K. 

The linearized dimensionless MHD equations with 
these assumptions are: 



dvy_ 
dt 



dby_ 
dt 



1 



Poix,z) 



Bo,j:{x, z)—^ + Boz{x, z)—^ 



-Vq' 



dz 



" ' v^^v. 



y 



(14) 



Box{x,z)-^ + Bo,{x,z)-^ 

dby „9, 

az 



(15) 



where densities, velocities, magnetic field, time and 
space coordinates are normalized to po (the plasma den- 
sity at dimensionless z = 6), Vao, Bq, t, and a (spicule 
radius), respectively. Also the gravity acceleration is 
normalized to a} jr. Second terms in the left hand side 
of Eqs. [13] and [15] present the effect of steady fiows. 
Eqs. [141 and [15] should be solved under following ini- 
tial and boundary conditions: 



(9) z;y(x, z,i = 0) = Vao exp 



2^ d ' 



by{x,z,t = 0) = 0, 



sin(fc2)e"/''^ 
(16) 



where d is the width of the initial packet. 

Vy{x = 0, z, t) — Vy{x = 2, Z, t) — 0, 
by(x = 0, z, t) ~ by{x = 2, z, t) — 0, 
Vy{x, z ~ 0,t) — Vy{x, z = 8, t) =: 0, 
by{x, z = 0, t) = by{x, z = 8, t) = 0. 



(17) 



Figure [T] shows the initial wave packet given by Eq. [16] 
for d = 0.5a (a is the spicule radius). The parameter k is 
chosen in such a way to have the standing Alfvcn wave. 



3 Numerical results and discussion 

To solve the coupled Eqs. [TU and [15] numerically, the 
finite difference and the Fourth-Order Runge-Kutta 
methods are used to take the space and time deriva- 
tives, respectively. We set the number of mesh- grid 
points as 256 x 256. In addition, the time step is cho- 
sen as 0.0005, and the system length in the x and z di- 
mensions (simulation box sizes) are set to be (0,2) and 



(0,8). The parameters in spicule environment are as fol- 
lows: a = 1000 km(spicule radius), d ~ 0.5a = 500 km 
(the width of Gaussian paeket), L = 8000 km (Spicule 
length), vq — 25 km s~^, rig = 11.5 x 10^^ m~^, 
Ba = 1.2 X 10-3 Tesla, To = 14 000 K, g = 272 m 
s-'^, R — 8300 m^s-^k"^ (universal gas constant), 
Vao = 75 km/s, fi = 0.6, r = 13 s, po = 1-9 x lO^^o kg 
m"^ po ^ 3.7 x 10^2 N m-^ ^o = 47r x lO^^^ Tesla m 
A"^, Zr = 6000 km (reference height), z^, = 200 km, 
zt = 2000 km, xo = 1, zq = 0.5, a = 2, H = 1000 km, 
ry = 10'^ m^ s~^, fcf, ~ 7r/8, an d fc = 7r/3 (dimensio nless 
wavenumber normalized to a) (JEbadi et al.ll2012bf ). 
Figure [2] shows perturbed velocity variations with re- 
spect to time in X = 1000 km, z = 1000 km; x ~ 
1000 km, z = 4000 km; and x = 1000 km, z = 7000 km, 
respectively. In Figure|31 perturbed magnetic field vari- 
ations are presented for x = 1000 km, z = 1000 km; 
X = 1000 km, z = 4000 km; and x = 1000 km, 
z = 7000 km, respectively. In these figures the per- 
turbed velocity and magnetic field are normalized to 
Vao a-nd Bq respectively, and it is obvious from the 
plots that there is a damping at the first stage of phase 
mixing. This behavior can be related to the presence 
of transition region between chromosphere and corona. 
At the first height {z = 1000 km), total amplitude of 
both velocity and magnetic field oscillations have val- 
ues near to the initial ones. As height increases, the 
perturbed velocity amplitude does increase in contrast 
to the behavior of perturbed magnetic field. Nonethe- 
less, exponentially damping behavior is obvious in both 
cases. This means that with an increase in height, 
amplitude of velocity oscillations is expanded due to 
significant decrease in density, which acts as iner- 
tia against oscillations. Similar results are observed 
by time-distance an alysis of solar spicule oscillations 
( Ebadi et al.ll2012al) . It is worth to note that the den- 
sity stratification influence on the magnetic field is neg- 
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ligible, which is in agreement with Solar Optical Tele - 
scope observations of solar spicules (jVerth et al.ll201l[) . 

Figures |4] and [5] illustrate the 3D plots of the per- 
turbed velocity and magnetic field with respect to x, z 
for t = IOt s, t = 50r s, and t = 80t s. In the presence 
of the transition region and stratification due to gravity, 
the damping process takes place in time than in space. 
It should be noted that due to the initial conditions, 
the damping time scale for the velocity field pattern 
is longer than the magnetic field one. In spite of the 
standing waves considered here, propagating waves are 
stable a nd dissipate after some periods due to phase 
mixing |Ebadi et al.ll2012U ). 

In Figure[Sl kinetic, magnetic and total energies nor- 
malized to the initial total energy, are presented from 
top to bottom respectively. Since spicules have short 
lifetimes and they are transient phenomena, we claim 
that in such circumstances th e phase mixing can occu r 
in time rather than in space ( De Moortel et al.lll999l) . 
Obtained damping times are in agreer nent with spicule 
lifetimes ( Zaaarashvili fc Erdelvill2009 ). Comparison of 
our results with the results deduced from the work of 



Fig. 1 (Color online) The initial wave packet for d = 0.5a 
is showed. 



Ebadi fc Hosseinpouij (|2013l ) shows that density gradi- 
ents because of transition region make the enhanced 
phase mixing in standing Alfven waves. As an exam- 
ple, the total energy amplitude damped to 10% of initial 
amplitude in 195 and 325 seconds, respectively. This 
means that the transition region between chromosphere 
and corona can accelerate damping rates of standing 
Alfven waves in solar spicules. 



4 Conclusion 

In spite of propagating Alfven waves, standing Alfven 
waves can also be found in open structures for example 
in spicules, and so may represent an efficient heating 
mechanism in the solar corona. In this paper, we as- 
sume that spicules are small scale structures in com- 
parison with coronal holes and other mega structures. 
In our model, we consider spicules with steady flow 
and sheared magnetic fleld. In addition, the transi- 
tion region between chromosphere and corona has also 
been considered. The considered medium is dense in 
its lower heights while becomes rare and rare as height 
increases (stratification due to gravity). Moreover, the 
density changes across the spicule axis which lead to 
different Alfven speeds (phase mixing). Our simula- 
tions show that phase mixing can occur in time rather 
than in space for standing Alfven waves. The per- 
turbed velocity amplitude increases as height is in- 
creased (JEbadi et al.ll2012a[ ). In contrast, the perturbed 
magnetic field amplitude decreases with height. Both 
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Fig. 2 (Color online) The perturbed velocity variations 
are showed with respect to time and x — 1000 km for three 
values oiz = 1000 km, z = 4000 km, and z = 7000 km from 
top to bottom. 




Fig. 3 (Color online) The perturbed magnetic field vari- 
ations are showed with the same coordinates as inferred in 
figure [2I 
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Fig. 4 (Color online) The 3D plots of the transversal com- 
ponent of the perturbed velocity with respect to 2:, z in 
t — lOr s, t = 50r s, and t = 80r s for kt = 7r/8. 



Fig. 5 (Color online) The same as in Fig. |4]for the per- 
turbed magnetic field. 
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Fig. 6 (Color online) Time variations of the normalized 
kinetic, magnetic and total energies for d = 0.5a. 



perturbed velocity and magnetic field amplitudes de- 
crease efficiently with time. Since spicules are short liv- 
ing and transient structures, the fast dissipation mech- 
anism is needed to deliver their energy to the corona. 
Our findings illustrate that the transition region can 
be accounted as a factor that can accelerate damping 
rates. We observed that kinetic, magnetic and total 
energy decrease with time exponentially. 

Acknowledgements This work has been supported 
financially by the Research Institute for Astronomy and 
Astrophysics of Maragha (RIAAM) , Maragha, Iran. 



References 

Alfven, H.: M.N.R.A.S. 107, 211 (1947) 

Aschwanden, M.J.: Physics of the solar corona. 

Beckers, J.M.: Sol. Phys. 3, 367 (1968) 

Cargill, P. J., Spicer, D. S., Zalesak, S. T.: Astrophys. J. 
488, 854 (1997) 

Del Zanna, L., Schaekens, E., Velli, M.: Astron. Astrophys. 
431, 1095 (2005) 

De Pontieu, B., Mcintosh, S.W., Carlsson, M., et al.: Sci- 
ence 318, 1574 (2007) 

De Moortel, I., Hood, A.W., Arber, T.D.: Astron. Astro- 
phys. 346, 641 (1999) 

Ebadi, H., Zaqarashvili, T.V., Zhelyazkov, I.: Astrophys. 
Space Sci. 337, 33 (2012a) 

Ebadi, H., Hosseinpour, M., Altafi-Mehrabani, H.: Astro- 
phys. Space Sci. 340, 9 (2012b) 

Ebadi, H., Hosseinpour, M.: Astrophys. Space Sci. 343, 11 
(2013) 

Erdelyi, R., Ballai, I.: Astronomische Nachrichten 328, 726 
(2007) 

Heyvaerts, J., Priest, E.R.: Astron. Astrophys. 117, 220 
(1983) 

HoUweg, J.V.: Astrophys. J. 181, 547 (1973) 

HoUweg, J.V.: J. Geophys. Res., 91, 4111 (1986) 

lonson, J. A.: Astrophys. J. 226, 650 (1978) 

Karami, K., Ebrahimi, Z.: Publ. Astron. Soc. Aust., 26, 
448 (2009) 

Kudoh, T., Shibata, K.: Astrophys. J. 514, 493 (1999) 

Kukhianidze, V., Zaqarashvili, T. V., Khutsishvili, E.: As- 
tron. Astrophys. 449, 35 (2006) 

McKenzie, J.F., Banaszkiewicz, M., Axford, W.I.: Astron. 
Astrophys. 303, 45 (1995) 

Matthaeus, W.H., Zank, G.P., Oughton, S., MuUan, D.J., 
Dmitruk, P.: Astrophys. J. Lett. 523, 93 (1999) 

Moriyasu, S., Kudoh, T., Yokoyama, T., Shibata, K.: As- 
trophys. J. 601, 107 (2004) 

Murawski, K., Zaqarashvili, T.V.: Astron. Astrophys. 519, 
A8 (2010) 

Okamoto, T.J., De Pontieau, B.: Astrophys. J. Letters 
736:L24, 6 (2011) 

Priest, E.R.; Solar magnetohydrodynamics. Reidel, Dor- 
drecht (1982) 

Ruderman, M.S., Berghnians, D., Goossens, M., Poedts, S.: 
Astron. Astrophys. 320, 305 (1997) 



Smith, P.D., Tsiklauri, D., Ruderman, M.S.: 2007, Astron. 

Astrophys. 475, 1111 (2007) 
Sterling, A.C. 2000, Sol. Phys. 196, 79 (2000) 
Verth, G., Goossens, M., He, J.S.: Astrophys. J. Lett. 733, 

15 (2011) 
Zaqarashvili, T.V., Erdelyi, R.: Space Sci. Rev. 149, 335 

(2009) 
Zaqarashvili, T. V., Khutsishvili, E., Kukhianidze, V., 

Ramishvih, G.: Astron. Astrophys. 474, 627 (2007) 



This manuscript was prepared with the AAS IATJ5X macros v5.2. 



